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The chemistry of molecular clouds has been studied for decades, with an increasingly general and sophisticated treatment 
of the reactions involved. Yet the treatment of turbulent diffusion has remained extremely sketchy, assuming simple Fickian 
diffusion with a scalar diffusivity D. However, turbulent flows similar to those in the interstellar medium are known to 
give rise to anomalous diffusion phenomena, more specifically superdiffusion (increase of the diffusivity with the spatial 
scales involved). This paper considers to what extent and in what sense superdiffusion modifies molecular abundances 
in interstellar clouds. For this first exploration of the subject we employ a very rough treatment of the chemistry and the 
effect of non-unifom cloud density on the diffusion equation is also treated in a simplified way. The results nevertheless 
clearly demonstrate that the effect of superdiffusion is quite significant, abundance values at a given radius being modified 
by order of unity factors. The sense and character of this influence is highly nontrivial. 
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1 Introduction 

The chemistry of molecular clouds is a complex system af- 
fected by different factors. The first models laid down the 
foundations of the chemical reaction network: following the 
initial steady state gas-phase chemistry model of Herbst & 
Klemperer (1973) several (pseudo-)time-dependent models, 
with fixed profiles of physical parameters were developed, 
considering more and more chemical reactions, while ne- 
glecting photodissociation (Leung et al. 1984; Herbst & Le- 
ung 1989; Millar & Herbst 1990). The agreement of the mo- 
dels with observed fractional abundances improved when 
the effect of turbulent diffusion was taken into account, re- 
sulting in smoother density distribution profiles for the more 
important species (Xie et al. 1995, hereafter XAL95). Sub- 
sequent models further refined the chemistry, considering 
details of the ion-neutral reaction scheme, and taking into 
account the gas-grains interaction, H2 and CO self-shield- 
ing, grain accretion effects and adsorption onto grains (Will- 
acy et al. 2002; Yate & Millar 2003). 

In contrast to the enormous efforts made in order to im- 
prove the representation of chemical processes, the treat- 
ment of turbulent diffusion has remained extremely sketchy, 
assuming simple Fickian diffusion with a scalar diffusiv- 
ity D. Under the conditions prevailing in interstellar space 
this is surely wrong. Interstellar turbulence is mainly driven 
by supernova shocks on spatial scales far exceeding that 
of the molecular clouds. From those scales, kinetic energy 
cascades down to the scale of clouds, cloud clumps and 
even cloud cores (Boldyrev 2002; Scalo & Elmegreen 2004; 
Ryan Joung et al. 2006). Indeed, in current thinking clouds 
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are hardly more than positive density fluctuations in this 
large-scale compressible turbulent flow. In such a turbulent 
flow, the increase of the rm.s. separation A of two tracer 
particles will only be due to eddies smaller than A, lead- 
ing to a scale-dependent eddy diffusivity. As a result, the in- 
crease of A with time will deviate from the Fickian square- 
root law: such non-Fickian or anomalous diffusive process- 
es commonly occur in various fields of physics (Avellaneda 
& Majda 1992), though astrophysical applications have so 
far been limited to solar physics (Petrovay 1999). As in the 
turbulent cascade the diffusivity increases for larger spatial 
scales, the transport of passive scalars by inertial-range tur- 
bulence belongs to the subclass of superdiffusive processes. 



The work presented here is a first exploration of the ef- 
fect of superdiffusion on the chemistry of molecular clouds. 
Our approach is described in Section 2. In order to get a 
first impression of the importance of superdiffusive effects, 
we neglect the intricacies of the chemical reaction network 
and treat their net effect as a simple parametric source term 
in the diffusion equation. As anomalous diffusivity cannot 
be described by any partial differential equation in phys- 
ical space, we are forced to work in Fourier space, with 
D = D{k) a function of the wavenumber k. To alleviate 
the ensuing difficulties, in Section 3.1 we first consider a 
case where the density of the main component is uniform 
throughout the cloud. A semianalytical solution is possible 
in this case; the results indicate that there are significant and 
nontrivial differences between the diffusive and superdiffu- 
sive cases. Section 3.2 then considers what changes are to 
be expected if the assumption of a uniform cloud is dropped. 
Section 4 concludes the paper. 
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Fig. 1 Number density as a function of the radius in a 
uniform density cloud for a tracer whose distribution scales 
as no ~ 1/r in the nondiffusive case (solid Une). Dashed 
line: normal diffusive case; dotted line: superdiffusive case. 
The parameters are: Tc = lO^yrs, and D = 4 x 10^^ (A); 
10^3 (B); D = 1.6 X lO^^cm^s-i (C). 

2 Formulation of the problem 

Consider a spherically symmetric molecular cloud, with r 
the radial coordinate. Let riHs {f) denote the number density 
of hydrogen, rii (r, t) the number density and /j = the 
fractional abundance of a tracer molecule i. On the basis of 
simple mixing-length arguments (cf. XAL95), if turbulence 
is characterized by a single dominant scale I and velocity v, 
the diffusion equation for the tracer reads 

dtni{r,t) = DV[Vni{r,t) - fi{r,t)VnuAr)] (1) 
+Si{r,t), 

where Si{r,t) is the source/sink term, representing the 
chemistry reaction scheme, and D = (Iv) is the diffusion 
coefficient. In previous work the fiducial values v ^ 1 km/s 
and Z ~ 0.1 pc were used, yielding diffusivities on the order 

of 1022-1023 cm2/s. 

As here we focus on the diffusive term, the chemical 
source term will be represented very roughly as a simple 
relaxation term 

no{r) - n{r,t) 



S{r,t) = 



(2) 



throughout this paper. Here no (r) is the diffusionless equi- 
librium solution, while Tc is a relaxation timescale on which 
the solution would converge to no(r), were there no diffu- 
sive effects. On the basis of the results presented in XAL95 
the value of Tc can be estimated to lie in the range 10^-10^ 
years. 

In a turbulent cascade regime where turbulence cannot 
be characterized by a single dominant scale anymore, equa- 
tion (1) does not hold. However, its Fourier transform, with 
the appropriate scale-dependent diffusivity value D{k) will 
still hold for each Fourier component of the passive scalar 
field. For the wavenumber dependence of the diffusivity we 
use the Kolmoeorov form D{k) = Da{k/ka)-^/^ through- 
out this paper, Dq being the value of the diffusivity at the 
wavenumber fco, corresponding to the fiducial scale I used 
in previous work. 
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Fig. 2 Same as in Fig. 1. but the parameters are: 
D = 1023cm2s-\ and Tc = 10^ (A); 7.6 x 10^ (B); 
1.2 x lOVs (C). 
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Fig. 3 Fractional abundance of O2 as a function of ra- 
dius in a uniform density cloud, with D = 1023cm2/s, 
Tc = 10^ yrs. Solid line: nondiffusive case, adopted from 
XAL95. Dashed line: normal diffusive case; dotted line: su- 
perdiffusive case. 

3 Results 

3.1 Uniform density cloud 

The Fourier transform of equation (1) takes a particularly 
simple form if the main constituent (i.e. H2) is distributed 
uniformly. Equation (1) then simplifies to 

dtn{r, t) = DV^n{r, t) + S{r, t) (3) 
The spatial Fourier transform of this is 

dth{k, t) = -D{k)k'^h{k, t) + S{k, t) (4) 
In the stationary case the solution is 

Transforming h{k, t) back to physical space we obtain the 
equilibrium distribution of the tracer. 

To carry out these transforms in the spherical geome- 
try at hand, we recall the theorem that an n dimensional 
Fourier transform can be replaced a one dimensional Han- 
kel transform, if the transformable function depends only 



© 2007 WILEY- VCH Verlag GmbH & Co. KGaA, Wemheim 



www.an-joumal.org 



Astroii. Nachr. / AN (2007) 



873 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 

Radius (10" cm) Radius (10" cm) 



Fig. 4 Same as in Fig. 3. but the inner and outer parts 
of the cloud are shown separately. The values of param- 
eters are: = lO^yrs, D = 4 x 10^^ (a), io23 (B) and 
D = 1.6 X lO^^cm^s-i (C). 

on r = x'f (Sneddon 1951). In three dimensions the 
transformation formulae are: 

n{k) = / r^n{r)Ji{kr)dr (6) 
Jo ^ 

n{r)= kih{k)Ji{kr)dk. (7) 
Jo ^ 

where Ji is the Bessel function of order h . These Hankel 
transforms can only be computed analytically for some spe- 
cial cases, such as a Gaussian profile for no(r) (cf. Mar- 
schalko 2006). For the more reaUstic cases presented below, 
equations (6) and (7) were evaluated numerically. 

In figures 1 and 2 we present analytical solutions assum- 
ing no(r) ~ r~^, a very rough analytical approximation for 
the actual profiles. It is apparent that there is a significant 
difference between the diffusive and superdiffusive case. In 
the superdiffusive case the concentration of the tracer in 
the cloud core is much more pronounced, as a consequence 
of the reduced diffusivity at small scales. Interestingly, an 
overall increase of the diffusivity D leads to an increase 
in the central number density of the tracer in the diffusive 
case, while it leads to a decrease in the superdiffusive case 
(though this trend reverses at large values of r). 

Next, instead of the generic l/r distribution we con- 
sider one particular tracer, the oxygen molecule. The form 
of no(r) for this molecule is taken from XAL95. Figure 3 
shows the fractional abundance of O2 as a function of the 
radius, still assuming a uniform H2 distribution. The diffu- 
sive and superdiffusive cases are again significantly differ- 
ent, but this distinction is better seen in the non-logarithmic 
plot (Figure 4). In contrast to the 1 /r case with its singular 
chemical source function, in the superdiffusive solutions the 
tracer is now less concentrated to the core than with Fickian 
diffusion. 

3.2 Non-uniform cloud 

The non-uniformity of the cloud was in fact implicitly 
already taken into account in the previous subsection, as 

www.an-journal.org 
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Fig. 5 Comparison of the O2 fractional abundance pro- 
files computed assuming normal (Fickian) diffusion with 
uniform (dotted) and non-uniform (dashed) cloud density. 
Case A:D = 10^3 cm^/s, = 7.6 x 10^ yrs. Case B: D = 
4 X 10^2 cm2/s, Tc = yrs. Case C: D = 10^3 cm^/s, 
Tc = 10^ yrs. SoUd: initial profile no(r). 

no(r) was taken from non-uniform cloud models. To make 
our approach more consistent, it would be necessary to take 
into account this non-uniformity also explicitly in equation 
(1). In this more general case the problem becomes math- 
ematically much more complicated. In order to get a first 
impression of the importance of the inhomogeneity, here we 
first compare the results obtained with and without explic- 
itly treating non-uniformity for the case of normal diffusion. 
For these calculations the density profile riHs (r) was taken 
from XAL95, as above. 

For the solution of equation (1) a time-relaxation method^ 
was used with a finite difference scheme first order accurate 
in time and second order accurate in space, on a grid of 
512 points evenly distributed between rin = 4 x 10^^ and 
''out = 4 X 10^^ cm. At the inner boundary of our domain 
the boundary condition was drn{r = rin,t) = 0, while 
at the outer boundary the number density was set to zero 
n(r = Tout 5 i) = 0. Starting from the same initial con- 
ditions as in the previous case, the solution was allowed 
to evolve in time. Comparing the solutions obtained in the 
uniform case to the exact stationary solution obtained from 
equation (5) we found that satisfactory agreement is reached 
within one characteristic time Tc. 

The results are plotted in figure 5. The effect of non- 
uniformity is clearly to increase the tracer abundance, espe- 
cially in the outer parts of the cloud. Nevertheless, vary- 
ing the parameters of the problem {Dq, Tc) the resulting 
effects are qualitatively similar and quantitatively compa- 
rable to those found in the uniform case. This encourages 
us to attempt to characterize the effect of superdiffusion in 
the nonuniform case by rescaUng the uniform superdiffu- 
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Fig. 6 Fractional abundance of O2 as a function of radius 
in a cloud of non-uniform density, with D = 10^^ cm^/s, 
Tc = 7.6 X 10^ yrs. Dashed line: normal diffusion. Dot- 
ted line: the expected abundance profile with superdiffusion, 
obtained by the rescaling method described in the text. 



sive solution for the fractional abundance (dotted line in fig- 
ure 3) by the ratio of the nonuniform and uniform solutions 
in the normal diffusive case. The resulting abundance curve 
is plotted in figure 6. Just as in the uniform case (figure 3), 
superdiffusion reduces the O2 abundance in the cloud core 
and increases it in the envelope. 



4 Conclusion 

In this paper we considered to what extent and in what sense 
superdiffusion modifies molecular abundances in interstel- 
lar clouds. The results clearly demonstrate that the effect is 
quite significant, abundance values at a given radius being 
modified by order of unity factors. The sense and character 
of this influence is highly nontrivial. 

In order to make a meaningful comparison with obser- 
vational data possible, this model needs to be developed fur- 
ther to include a realistic treatment of the chemical reactions 
and to include a more rigorous calculation of the effect of 
non-uniformity of the cloud in the diffusion equation. Work 
in this direction is in progress. 
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